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Abstract 


We calculate the 0(a s ) corrections to the charged Higgs mediated inclusive 
semi-tauonic B decay. By working in the Landau gauge, we demonstrate 
how to obtain the charged Higgs contributions (both direct and interference 
terms) from the known QCD corrections to t —» bH + . Combining our re¬ 
sults with 0(l/m?) corrections to the spectator model and comparing the 
theoretical prediction with the recent experimental measurements, we find a 
model-independent 2cr upper bound on the ratio of Higgs vacuum expectation 
values, tan/3 < 0.52 (m#/l GeV). 
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I. INTRODUCTION 


The B —> X t v decay is sensitive to extensions of the standard model (SM) Higgs sector. 
In particular, for two Higgs doublet models (2HDM) |l|, the experimental measurement 
of this decay provides the strongest upper bound on tan [3/mH, when the ratio of Higgs 
vacuum expectation values, tan /3, is large [^j-||]. Since the charged Higgs boson contributes 
at tree level, its contribution cannot be canceled by other new particles in the theory. 
Consequently, this bound holds in all SU(2)xU(l) clectroweak models with a charged Higgs 
boson, independently of other details of the model. It is therefore important to calculate 
the charged Higgs contribution with high accuracy. 

The leading order calculation is based on tree diagrams of the spectator model. Correc¬ 
tions to spectator model results are given by the 1 /mb expansion in heavy quark effective 
theory (HQET) ||. The 0(l/mf) corrections have been computed in the case of massless 
leptons p-|8|, and for the tau-lepton in both the SM P~ m| and the 2HDM P|,|T2f|. The 
one-loop QCD corrections to the tree-level amplitudes of the spectator model are equally 
important. In the SM, the 0(a s ) corrections have been computed for both massless lep¬ 
tons mm and the tau-lepton |T6|-|T8|. However, no such calculations are available for the 
2HDM. Their absence is a major source of uncertainty in the branching ratio predictions for 
large tan (3 || . 

The 0(a s ) corrections for t —> bW [[19] and t —> bH |2(J are known. The QCD corrections 
for the B —> X t v decay can be directly obtained from these results. This method was 
applied in ]17| to the SM calculation. In this paper, we show that by working in the Landau 
gauge, the method of |L7| can be extended to include the charged Higgs contributions. 

This paper is organized as follows. In section II we calculate the 0(a s ) corrections to 
the B —> X t v decay in the 2HDM (with model-II Higgs-fermion couplings [|lj). In section 
III we combine this new result with the known 0(1/ml) corrections and the quark mass 
relations of HQET, and get an accurate prediction for the branching ratio as function of 
tan/5 and the charged Higgs mass, m#. Section IV contains our summary. Improvements 
to the QCD computation through renormalization group resummation are discussed in the 
Appendix. 


II. 0(a t ) CORRECTIONS TO B=*Xtv 


In |7 


the QCD corrections to B —> X r v in the SM were obtained using the known 
QCD corrections of t —> bW and t —> bH decays, and without calculating any new Feynman 
diagrams. The idea is to decompose the n body phase space (where n=3 or 4) into a product 
of n — 1 body phase space and two-body phase space, integrated over the invariant squared 
mass of the tv pair 


21 


dPS(B ^Xtv) 


dq 2 dPS(b -> XW*) dPS(W* -> tv), 


( 2 , 1 ) 


where + p„ is the lepton pair momentum, X is either c or c+gluon, and dPS is 

the appropriate phase space differential element. The SM rate can be decomposed into two 
terms [|I7|: one corresponds to transverse W exchange and the other is equivalent to an 
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(un-physical) scalar exchange. The latter can be identihed with the exchange of a massless 
Goldstone boson in the Landau gauge. When generalizing the analysis of [|17]] to include 
a charged Higgs boson, there are two additional terms: one that arises purely from the 
charged Higgs exchange and the other that arises from the interference between the W and 
the charged Higgs diagrams. The method proposed in |17|] is now easily extended to this 


case by explicitly working in the Landau gauge. Then, there are three diagrams contributing 
to B —> X tv corresponding to VL-exchange, Goldstone boson exchange, and charged Higgs 
boson exchange. It is easy to show that the rate can be decomposed into the sum of two 
incoherent decays: 


M = A4\v + -M-g + -M-h 


\M\ 2 = \M w \ 2 + \M g + M 


H\ 


( 2 , 2 ) 


The proof of ( |2.2|) is straightforward. The lepton tensor in interference terms between a 
scalar (a Higgs or a Goldstone boson) and a vector boson ( W ) carries one Lorentz index, 
L fi . After integration over the tv phase space, L M must be proportional to the only Lorentz 
vector in the problem, q M (the lepton pair momentum). Contracting L M with the Landau 
gauge W propagator, T afl ~ g afl q 2 - q Q q we get 

L^T ail oc q^ig^q 2 - q a q^) = 0 . (2.3) 


This proves (ID- Note that this proof holds for arbitrary A" in Q2.1|) , and is therefore 
applicable to the QCD corrected rate to all orders in a s (at lowest order in the electroweak 
coupling expansion). 

Let us define 
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(2.4) 


Due to (|2.2|) we can write schematically (in the Landau gauge) 


“(!-Vp ) 2 dT 
dt 


dt , 


' Pr 


T(B -*• Xtv) = 

( _ = H «P(p ^ cW *) T afl Tp u ( W* 


(2.5) 


tv) + H(b -► cS*)L(S* -> tv) , 


where H and L refer to the terms that arise from analyzing the hadronic side and the 
leptonic side of the decay process, respectively. As above, T afM is the Landau gauge W 
propagator. The two terms on the right hand side of the above equation correspond to the 
sum over the two incoherent decays in (|2.2| ). In particular, S stands for an effective scalar 
that we construct to give the contributions of both the Goldstone boson and the physical 
charged Higgs. Moreover, it is convenient to absorb the Goldstone boson and charged Higgs 
propagator factors (—1 /q 2 for the Goldstone boson and 1 /m? H for the physical Higgs boson 
in the approximation that q 2 -C m 2 H ) directly into the definitions of the effective couplings 
of S to the leptons. As a result, the scalar propagator factors do not explicitly appear in 


( |2.5|) . The couplings of S to quarks are then given by 


igm b V cb 

2\/2m w 


a+ 675), 


( 2 . 6 ) 
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where a = a G + an, and b = b G + bn, with |jj 


a G ~ “(I — %fP)l t) G — —(1 + y/p), 

a H = tan [3 + yfp cot/3, b H = tan/? — yfp cot (3. (2.7) 


The couplings of S to leptons (with the propagator factors included as indicated above) are 
given by 


W 


2\f 7 2rnwrrib 

where a £ = a G + a e H , and b e = b G + b e H , with 

\/~Pt 


{a 1 + 6 ^ 75 )- 


a G ~ 


U - 

°G ~ 


y/Pr 


t 


a H ~ 


yfp^ tan f3 


,0 _ -y/fr tan/3 
b H — ~ 


t n z 

Integrating over the leptonic two-body phase space yields L [see (p.5| )j 


m 


1 - y ) tpr 


1 2 tan B tan 2 3 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


We denote T a fe c = T(a —> be). The expression for V t bH can be found in [12^1. For arbitrary 
quark flavors and scalar Yukawa couplings, Eqs. (7) and (8) in [^0| can be parameterized as 
rQg5(m 2 /r?7.Q, Ci, C2, C3). The relevant combinations of couplings are 


ci = a 2 + b 2 = 2(1 + p) + 2(tan 2 /3 + p cot 2 (5) — 4(tan (3 — pcotfi), 

C 2 = a 2 — b 2 = Ay/p [—1 + tan /3 cot (3 + (tan /3 — cot (3 )\, 

c 3 = ab = (1 — p) + (tan 2 (3 — pcot 2 (3) — 2 (tan /3 + pcot (3 ). (2-11) 


Then, the contribution to the B —> X t v decay rate from the effective scalar S is obtained 
by taking S to be off-shell and integrating over its momentum: 


rV-VP) 

Ts(B^Xtu)^ / FbcsiPi t', Ci, c 2 , c 3 ) L(t) dt , 

J p T 

where L(t) is given in ( p. 10| ) and Ci, c 2 and c 3 are given in (|2.11|) . 


( 2 . 12 ) 


III. NUMERICAL ANALYSIS 


In our numerical analysis we set cot f3 = 0, as we are always interested in large tan/3, 
and we neglect the tiny b —> u transitions. Since Tq q s is linear in the three relevant terms 
in (2. 11|) we can identify the Higgs and the interference terms in the total rate 
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dt. 


(3.1) 


where T# (T/) is the Higgs (interference) term in E The tree level decay rate can be obtained 
by inserting the appropriate tree-level expression for T bc s info ( P • 1|) • We have performed the 
integrals analytically and checked that we recover the known results [[|.fjj. Turning to the 
0 (a s ) computation, we insert the appropriate one-loop expressions for F bcS deduced from 


The integrals in (|3.1|) are then evaluated numerically. 

Our results are not very sensitive to small variations in the choice of the c and b pole 
masses. Below we quote results for representative values of the pole masses m c = 1.4 GeV 
and m b = 4.8 GeV. Decomposing the total rate into the tree-level and the one-loop parts, 


T = T° + a^T 1 , 


(3.2) 


it is convenient to quote results for T 1 /T° (note that factors of tan/5 cancel in the ratio). 
We get 

pi pi 

iif - _n 7^3 iI 

r° H po 


p 1 pi 

tXL = _o 450 t± 

1 W 1 e 


= -0.487. 

(3.3) 

on the calculations of 

16 , and in 17 : 

= -0.545, 

(3.4) 

Xtv (B -> Xev). 

Defining p = T/T 0 , 

7 ) = 0.115], we get 


0.89, r) H — 0.83. 

(3.5) 


The main uncertainties in the calculation are from the unknown 0(a 2 ) corrections. To 
account for these, we vary a s in the range 0.20 < a s (p) < 0.36, corresponding to m b /3 < 
p, < m b and 0.110 < a s (mz) < 0.125. This leads to the ranges 


1.02 < -BL < 1.04, 1.01 < — < 1.03, 0.96 > — > 0.91. (3.6) 

be Ve Ve 

Note that these ratios are rather insensitive to the exact choice of a s . This completes the 
analysis of the 0 (a s ) corrections. 

To arrive at our final prediction, one must also include the 0(l/m 2 ) corrections to the 
spectator model. In the Heavy Quark Expansion, we use the HQET mass relation 


T Ai + 3 A 2 

m B = m b + A--- 

2 m b 


+ ... , 


rrio — rn c + A — 


Ai + 3 A 2 

2 m c + 


(3.7) 


and the numerical values as in [[| 
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0.4 < A < 0.6 GeV, 0 < -Ai < 0.3 GeV 2 , 0.11 < \ 2 {m b ) < 0.13 GeV 2 . (3.8) 


with our 0(a s ) calculation, and using the range 
, we get the theoretical prediction 


Combining the 0(1/ml) calculation of 
of the various parameters as given in Eqs. ( |3.6| ) and 
for the branching ratio (normalized to the electron channel) as a function of r = tan (5/mu 
given in Fig. 1. Comparing it with Fig. 1(a) of [£|] we learn that, as anticipated in |4|], the 
errors in the plot are reduced by about a factor of two for large r. The residual errors arise 
almost entirely from the unknown pole masses, namely the values of A and Ai. 

Combining the recent measurements of the branching ratio of B —> X r v from ALEPH 
and L3 and using the world average for B —> X £ v (^ = e, n) P4fl , 

(3.9) 


BR (B -> X tv) = 2.69 ± 0.44 1 


BR (B -> X i v) = 10.43 ± 0.24' 


we obtain the 2cr upper bound 


r < 0.52 GeV 


-i 


(3.10) 


The ler upper bound is r < 0.49 GeV -1 (as compared to r < 0.51 GeV -1 in 0). The improve¬ 
ment over the bounds obtained in 0 originates from both the experimental improvement 
and from our 0(a s ) calculation. 

Our analysis holds with minor modifications for general multi-Higgs doublet models with 
natural flavor conservation (for a recent analysis, see [^5|). In such models, instead of the 


single parameter tan /?, three complex coupling constants determine the Yukawa interactions 
of the lightest charged scalar (assuming that the heavier charged scalars effectively decouple 
from the fermions). The parameters X, Y and Z describe the couplings to down-type 
quarks, up-type quarks and charged leptons, respectively. The bounds on \J\XZ*\/mu and 
\J\Re(XZ*) \/mu are the same as that on r. The 2 a upper bound on the imaginary part of 
XZ* is \J |Im (XZ*) |/ m H < 0.44GeV" 1 (< 0.41 GeV” 1 at la). 


IV. SUMMARY 

We have calculated the 0(a s ) corrections to B —> X r v decay in the 2HDM. In addition 
to the known SM result, we have computed the contributions arising from the exchange of 
a charged Higgs boson; both the pure Higgs term and the W + ~H + interference term have 
been included. Working in the Landau gauge, the calculation is carried out by making use of 
previous calculations of the QCD corrections to t —> bH. Our results for the numerical value 
of the QCD corrections are given in Eq. (|3.5|) . We combined our calculation with the known 
0(l/ml) corrections and the quarks mass relations of HQET, and obtained an accurate 
prediction for the branching ratio as function of tan shown in Fig. 1. The remaining 

uncertainties in the prediction are almost entirely from the uncertainties in the pole masses. 

The virtual effects of the charged Higgs boson can be important in other processes. 
However, in most cases, the charged Higgs exchange only first enters at the one-loop level. 
In models with an extended Higgs sector, there is often additional new physics beyond the 
Standard Model that can also contribute to the processes under consideration. In this case, 
there can be no clean bound on Higgs sector parameters, since the one-loop charged Higgs 
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effects can be partially canceled by the effects of other particles in the theory. For example, 
in the minimal supersymmetric model, the bounds on the charged Higgs parameters from 
the one-loop charged Higgs contributions to b —> sy |2(| and Z —> bb |^7| decays can be 


avoided if partially canceled by loop diagrams with intermediate superparticles. In contrast, 
the charged Higgs boson contribution to B —> A" r v studied in this paper occurs at tree 
level; therefore, its effect cannot be canceled by other new particles in the theory. 
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APPENDIX A: RENORMALIZATION GROUP IMPROVEMENT 


The QCD corrections to T(t —> bH) were calculated for = 0 in and extended 

Hence, we use the results of |20||[*| . All calculations were performed 


for finite in 20,30 


in the on-shell mass renormalization scheme, in which the quark pole masses appear in the 
formulae. These results can be improved by using the technology of the renormalization 
group equations (RGE) to resum the leading logarithms to all orders in a s . In particular, 
note that for small e = m 2 /m 2 , a,,loge is typically of 0(1). Consequently, all leading 
logarithmic terms of the form ea™log n e should be summed j31|. This can be accomplished 
by noting that one of the effects of the one-loop QCD corrections is the replacement of 
the pole mass mb with the one-loop running mass mb(mt) |32|. By replacing the one-loop 
mb(m t ) with the fully RGE-integrated running mass, m EGE (m t ), one successfully resums all 
a s log e terms. 

The one-loop QCD corrections to T(t —> bH) were obtained in [|20|, but the RGE- 
improved rate was not given. To carry out the RGE-improvement, first expand the results 
of [Q in powers of e. Formally, the full QCD-corrected rate can be written as 




n __ m-pm 

€ a c L „ 


(Al) 


which defines T™. In this case, only terms up to 0(a s ) are known. It is sufficient to keep 
terms up to 0(e) since e ~ 10 -3 is a small number, and only the 0(e) term can be enhanced 
by the possibly large tan 2 /? coefficient. Thus, we write V = rQ + a s rQ + e (T) 1 + a s r(). Using 
the expressions given in [[2(J|, we find 


*The results of [ 301 differ from those of [|(J. However, the latter results appear more reliable; see 
the discussion of [2C] for further clarifications. 
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r = r° + e r° 


+ a s r /m> 


(A 2 ) 



4«s 

71 


loge 


where Ty- in contains no loge terms, and 
G F ml\V tb \ 2 


r° — 
1 0 — 


r? = 


8V2n 

G F m 3 \V tb \ 2 

2y/2: 


(l-x)cot (3, 


71 


(1-x) 


1 + ^(1 - x) tan 2 fj 


~ | cot 2 p }> , 


(A3) 


where x = rrr H /rrq. Using 


m 2 b (m t ) = ml(m b ) ^1 + loge + 0 (ct 2 )^ , 


(A4) 


we see that (|A2| ) is equivalent at 0(a s ) to replacing m b (m b ) with m b {m t ) in To + eT) 1 . This is 
a general result [^]. Then, to sum up the ect”log n e terms, we replace the 0(a s ) expression 
for m b (rrit) with m^ GE (? 7 it) [|31||: 


2a, 


m b (m t ) = m b (m b ) ( 1 + loge 




m EGE {m t ) = m b (m b ) 


ois{jn b ) 

a s (m t ) 


- —2"/o/Po 


(AS) 


where 70 = 2 and Pq = 11 — 2/3 rif. 

The decay rate depends on m b through both the Yukawa coupling and the phase space 
integrals. Far from threshold, the dependence of the phase space integrals on m b is weak; 
therefore, the RGE-improvement has no numerical significance there. In contrast, the RGE- 
improvement via the Yukawa couplings is important provided that tan P is large (it is only 
for large tan P that the rate is sensitive to m b ). To see the effect, we define 7 = T/T° 
prior to RGE-improvement (see section 111), where in the notation above, T° = 

The results of the RGE-improvement are sensitive to the procedure that we take. We can 
use the RGE-improved m b (m t ) in T° only (we call the result t]q GE ) or in both T° and T| m 
p]i GE ). The difference is 0(a 2 ), beyond the accuracy of the calculation, but in practice it 
can turn out to be important. We take m b = 4.8 GeV, m t = 175 GeV, rn h = 100 GeV 
and a s {mt) = 0.1. As a first example, take tan P = 1 . Then 7 ~ rjp GE ~ ^RGE ^ 0.92, 
with difference of O(10~ 4 ). We see that the weak dependence of the rate on m b makes the 
RGE-improvement unimportant. As a second example, take tan P = 100: 


r) ~ 0.46, i] EGE ~ 0.50, r] EGE ~ 0.46. (A 6 ) 

The difference between r/ EGE and r)[ iGE is a measure of the size of the neglected 0(a 2 ) terms. 
Given that 77 ~ r] EGE , it follows that for very large tan/?, the RGE-resummed result does 
not provide a practical improvement over the original one-loop calculation. 

The situation is different for the B —> A 'tv decay in three important ways: 

(i) The RGE improvement is related to the Yukawa coupling of the light quark. In 

B X'tv, the Yukawa coupling of the charm can be safely neglected. 

(ii) The RGE improvement is important to sum large logs. However, log (m b /m c ) 0 ( 1 ) 

is not large. 
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(iii) The RGE improvement is a correction to the 0(a s ) corrections, but as noted below 
Q3.6|) , the latter almost cancel in the relevant ratio of rates that we consider. 

Therefore, while theoretically important, the RGE improved result has no numerical 
significance for the B —> X r v calculation. 
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FIGURES 
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FIG. 1. r(R — ► X t v) /T(B —* X € is') as a function of r = tan fi/mu • The shaded area 
between the solid lines is our result. The dash-dotted lines give the experimental 1 a bounds. 
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